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We study the effect of quenched random field disorder on a driven elastic interface close to the de- 
pinning transition at the upper critical dimension dc = 4 using the functional renormalization group. 
We have found that the displacement correlation function behaves with distance x as (In a;Ao)^'''^ 
for large x. Slightly above the depinning transition the force-velocity characteristics is described 
by the equation v ~ / |In/|'^^^, while the correlation length behaves as Lv ~ f^^^'^ |ln/|^''^, where 
/ = F/Fc — 1 is the reduced driving force. 
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The driven viscous motion of an elastic interface in a 
medium with randomly distributed pinning disorder has 
attracted considerable theoretical interest during the last 
decade and is in a state of rapid development. The rea- 
son is that, on one hand elastic interfaces in a disordered 
medium exhibit the rich behavior of glassy systems and 
on the other hand it can serve as a model for many exper- 
imental systems, such as domain walls in magnetically or 
structurally ordered systems with impurities and inter- 
faces between immiscible fluids in porous media. Other 
closely related problems are the motion of a vortex line 
in an impure superconductor and the motion of a dis- 
location line in a solid 0Q|- In all these systems the 
basic physical ingredients are identical: the elastic forces 
tend to keep the interface flat, whereas the disorder lo- 
cally promote the wandering. In the dynamics, this in- 
terplay between quenched disorder and elasticity leads 
to the complicated response of the interface to an exter- 
nally applied force. At zero temperature, a driving force 
F exceeding certain threshold value Fc is required to set 
the elastic interface in steady motion. The depinning 
transition can be considered as a nonequilibrium dynam- 
ical critical phenomenon where a system becomes ex- 
tremely sensitive to small perturbation. Recently, sig- 
nificant progress has been made in understanding the 
depinning transition 0,0 (for recent studies see jaQ). 
It has been shown that the functional renormalization 
group (FRG) gives an adequate description of the criti- 
cal behavior at the depinning transition if one presumes 
to consider a singular renormalized random force correla- 
tor. The scaling analysis shows that the disorder effects 
dominate over elasticity in dimensions c? < 4, and there- 
fore dc = 4 is the upper critical dimension of the problem. 
Below dc the interface undergoes the depinning transition 
at a critical driving force Fc and slightly above the criti- 
cal force Fc the average velocity of the interface behaves 
as 

Y^{F-Fcf, F>Fc, (1) 

where 6 being the critical exponent. The roughness ex- 
ponent characterizing the width w of the wandering in- 
terface at the depinning transition is defined by 

w ~ (2) 



where L is the linear size of the interface. The FRG anal- 
ysis carried out in Ref. p| enabled one to compute the 
critical force Fc, and the critical exponents 9 and C to or- 
der e = 4 — d. In the case of random force disorder it was 
found that = l-e/9 + 0{e'^) and C = e/3 HQ. It was 
suggested in that at the depinning transition the crit- 
ical exponents for random bond and random field (RF) 
disorder are the same. Initially, the result for the rough- 
ness exponent was expected to be exact for all ci = 1,2,3 
[3, but more recently, the non-zero two-loop correction 
to C lias been found u\- 

The results of FRG analysis have been checked by in- 
tensive numerical studies using both direct simulation 9- 
and different cellular automata models 13-15] , which 
are believed to belong to the same universality class. The 
computed values of critical exponents are in a good agree- 
ment with the predictions of FRG at least for d = 1,2,3. 
In the numerical works Il6j| the depinning transition 
was studied at the upper critical dimension. However, 
to our knowledge, no explicit consideration of the depin- 
ning transition at the upper critical dimension dc = 4 is 
available so far. Another motivation to consider the de- 
pinning transition at the upper critical dimension is that 
some experimental elastic systems, for example systems 
with dispersive elastic constant such as moving geologi- 
cal faults arising from earthquakes 0, or systems with 
long-range Coulomb interaction, have the upper critical 
dimension dc = 3 or 2. One expect that these systems 
may show a behavior which is similar to the behavior of 
a simple model at dc = 4 ^3 ■ 

It is well-known that at the upper critical dimen- 
sion the power laws modify to logarithmic corrections. 
While at the upper critical dimension the one-loop RG 
consists in summing the main logarithms, the two- loop 
RG takes into account the subdominant logarithms |19| . 
Due to the fact that close to the depinning transition 
the main logarithms are leading, the results of the one- 
loop RG treatment are expected to be exact at the upper 
critical dimension for F ^ Fe- 
in this Brief Report we consider the motion of an elas- 
tic interface in a disordered medium and our main pur- 
pose is to describe the critical dynamics near depinning 
threshold for d = dc by using FRG method to one-loop 
order. The motion of a d-dimensional interface obeys the 



2 



equation 



A 



dz{x, t) 
dt 



jy^z + F + g{x,z) 



(3) 



where A is the friction coefficient (or the inverse mobil- 
ity), 7 is the stiffness constant, and F is the driving force 
density. The quenched random force g{x, z) is assumed 
to be Gaussian distributed with the zero mean and the 
correlator 



{g{x, z)g(x' , z')) = S^"^^ {x ~ x')A{z — z'). 



(4) 



To make this model well-defined one has to introduce the 
cutoff A|^^ in the 5'^{x) function at scales of order the im- 
purity separation or other microscopic scales. We restrict 
our consideration to the case of random field disorder 
when the correlator A(z) = A(— z) is a monotonically 
decreasing function of z for z > and decays rapidly to 
zero over a finite distance. 

In Ref. Q the RG analysis of the model Q was 
carried out using the technique of path integrals in the 
one-loop approximation. After integrating out fluctua- 
tions in the momentum shell A; < |fc| < Aq, the following 
RG flow equations have been obtained 



d In A Kd f , f , ;/ ,~ X 

^ = "Tat y ^^^^ 



oo 



(5) 
(6) 



where v is the average velocity of the interface. A; — 
Aoe"', V = Av/(7A2) and K'^ = 2'^-^Ti'^''^T{dl2). Due 
to the tilt symmetry the stiffness constant 7 does not 
renormalize. In the limit v — > the disorder correlator 
A(z) renormalizes as 



dA(z) 
dl 



72 Af dz"^ 



A2(z) - A(z)A(0) 



(7) 



The RG equations Q and (O are the basis for computa- 
tion of the force- velocity characteristics in the vicinity of 
the depinning transition. In following we analyze Eqs. (0) 
and 10 at c? = 4, i.e. for e = 0. Before considering the 
general solution of the RG equations, we will analyze the 
flow equation for A"(0) 



dl 



(8) 



From Eq. ((SJ it follows that as in the case d < dc, the sec- 
ond derivative of the disorder correlator at origin A"(0) 
diverges at the finite length Ic = —j^/{3K4Aq{0)) for any 
initial condition Aq(0) < 0. Thus one obtains the Larkin 



length Lc = Aq e = at the upper critical dimension. The 
divergence of the curvature of A(z) implies the genera- 
tion of a cusp singularity: A(z) becomes non-analytical 
at the origin and acquires for I > lc & non-zero derivative 



A'(0"'") < 0. It was shown in Ref. that the cusp gener- 
ated during the renormalization determines the threshold 
force of the depinning transition. Therefore we expect 
that even at the upper critical dimension the interface is 
pinned for small enough driving force. 

We now will consider the depinning transition at the 
upper critical dimension. Although Eq. ||SJ| does not 
have a sense beyond the Larkin scale, nevertheless, we 
can still use the flow equation for the renormalized 
correlator. In contrast to d < dc where the critical be- 
havior at the depinning transition is obtained from the 
fixed-point solution of Eq. corresponding to the con- 
dition dA*{z)/dl = 0, the solution of Eq. describ- 
ing the behavior at the depinning transition at the up- 
per critical dimension depends explicitly on L To find 
the scaling form of the function Ai{z) at c? = 4 we 
look for an automodel solution of Eq. in the form 
Ai{z) = K^^j'^(t){l)p{ip{l)z). Note that the latter re- 
fiects the scaling behavior at the depinning transition. 
Substituting this scaling ansatz into Eq. we obtain 
the simultaneous equations for 4'{l), ip{l), and p{z) 

(t>'{i) = -4>HiW{i), ^'(0 = -0(O^'(O, (9) 

{zp{z))' ={p\z)/2-p[z)p{0))" . (10) 

Eqs. imply that 'iIj{1)/4>{1) = a is a constant which 
will be determined below. This condition allows us to 
find (j){l) = (Sa'^l)-^/^ and iP{l) = {3l/a)-^/^, so that the 
automodel solution of Eq. results in 



Ai{z) = K^^-f^{3a^l)-^/^p{z{3l/a)-^/^). 



(11) 



Eq. (fTT|l is the pendant of the fixed-point solution of the 
disorder correlator at d < 4. One should bear in mind 
that the ERG equation in this case {d = dc) gives the 
exact large-scale behavior while for d < dc one must rely 
on the e-expansion. The solution of Eq. (|10|l with the 
initial condition p{Q) = 1, which formally coincides with 
the equation for the fixed-point disorder correlator at e = 
3 can be written as 



p{z) - 1 -lnp(z) = z^/2. 



(12) 



Note that p{z) has a cusp at origin so that its behavior 
near z = is given by p{z) = 1 — |z| -I- ^z'^ + ... The 
constant a must be defined from the initial condition for 
the disorder correlator. Indeed, the flow equation 
for the disorder correlator implies that in the case of RF 
disorder the RF strength c — f^°° Aidz is conserved to 
one-loop order pl (it was shown in 0, Q that the above 
integral is not conserved in the two-loop order FRG), i.e. 
it does not depend on I. Therefore the constant a in the 
ansatz (|11|) is determined by the strength c of the bare 
disorder correlator as a « 1.55i4r^^7^c~^, where we have 
used the integral p{z)dz ~ 1.55. To higher orders 
of FRG the non-universal constant a is determined by 
higher moments of the bare disorder correlator. For the 
bare disorder correlator being a smooth function, the RG 
flow generates as in the case d < dc singularities on the 
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scale Ic, which result in the cusp of the running disorder 
correlator, and therefore in the existence of the threshold 
force Fc > 0. Because the automodel solution has 
the cusp on all scales, one should use it only beyond the 
Larkin scale. The full solution of the flow equation lO, 
of course, depends on the initial condition, nevertheless, 
the latter is expected to approach the solution Ullfl in 
the limit Z — > oo. From Eq. Hll|) one can immediately 
derive the scaling relations for the first derivatives of the 
running disorder correlator 



MO) 

a;(o+) 
Ar(o) 



2N|-1/3 



(13) 
(14) 
(15) 



Using the above results, we will now calculate the dis- 



placement correlation function B{q) 



[ZnZ- 



j) that de- 



scribes the roughness of the interface at the upper critical 
dimension. Simple scaling analysis shows that the corre- 
lation function satisfies the following fiow equation 



B{q)^e^'Biqe';Ai,Fi] 



(16) 



In order to extract the behavior for long- wavelength cor- 
relations at the depinning transition, F — Fc, we put 
qe' = Aq and expand Eq. (|16|l in powers of A. After 
some algebra this yields 



B{q) 



MO) 



1 



(3a2)V3 q4 [InAo/g]^/^' 



(17) 



where in order to obtain the final expression on the right- 
hand side of l|17|) we have used Eq. (|13|l . In a direct anal- 
ogy with the case d < 4 Eq. (|17|l holds simultaneously in 
the equilibrium and at the depinning threshold at least 
within the one- loop approximation [a, • The Fourier 
transform of Eq. (|17(l results in the following real-space 
displacement correlation function for large distance x 



B{x) - (lnxAo)^/^ 



(18) 



To obtain the force-velocity characteristics we have to 
integrate the flow equations in the vicinity of F^ i.e. in 
the limit of small v. Substituting Eq. (|15|l into the flow 
equation for the friction coefficient and Eq. H14|) into 
the flow equation for the driving force © we obtain 



din A, 



dl 
dFi 
'df 



2 

7^0 ,-2/3 -2i 



(9« 



1/3 



(19) 
(20) 



Eqs. H19(l and (|2()|1 describe the renormalization of A; and 
Fi to one-loop order beyond the Larkin scale at the up- 
per critical dimension. In contrast to Eqs. (|19|l and H20() 
below the Larkin scale the friction coefficient A/ in- 
creases under renormalization in accordance with Eq. jSJl, 
while the driving force is essentially not renormalized 



{dFi/dl ~ 0). Integration of Eqs. (jTH and over I 
starting from 1^ yields the following scaling relations for 
the friction coefficient and the driving force 

A, = Ao(^//c)-'/^ (21) 
Fi-F - -F, ~ -0.2A2(7c)1/3, (22) 

where Ao is the friction coefficient on the scale Ic- In 
order to obtain the renormalized friction coefficient one 
should express I in (|21|l through the correlation length 
Lv according to — Ag^^e'"'. In d < dc dimensions the 
relation between and v reads, ^ y-'^/i^-C) ^ where 
z = 2 — (e — 0/3 is the dynamic exponent relating the 
time and the length scale .4»5j . Let us now derive the 
relation between and v at the upper critical dimen- 
sion. To do this we need first the relation between the 
time scale and the space scale, which is derived by us- 
ing the relation t ^ XiL'^/'j with A; given by Eq. H21|) 
as t ^ L^(lnL)~^/^. Following '^l the correlation length 
Lv can be derived by equating the systematic drift of the 
interface and the height fluctuation (or equivalently from 
the velocity-velocity correlation function at equal times 
6]), which is the square root of Eq. ((TH|l . vt ^ (InL)^/''. 
Combining the latter with the above relation between t 
and L gives the correlation length at the upper critical 
dimension as 



.-1/2 



|lnv 



5/18 



(23) 



The use of the relation, Fi^ = vA;^ , which is obtained by 
stopping the renormalization at 1^, where Fi^ is given 
by H22I) . A/^ is given by (|21|l with 1^ = InLvAo, and 
Lv v^i/^ (the logarithmic correction to Lv in Eq. H23|) 
results in higher order terms in the force-velocity char- 
acteristics) gives the implicit form of the force-velocity 
characteristics in the vicinity of the depinning transition 
at the upper critical dimension as 



F-Fc 



|ln V 



,2/9 ■ 



(24) 



Within the one-loop consideration the interface velocity 
V under the logarithm of Eq. (|24|l can be replaced by 
F - Fc. Substituting v from Eq. into Eq. we 
express the correlation length as function of the driving 
force as Lv ~ /~^^^ |ln/|^^^, where the reduced driving 
force f = F I Fc — 1 is introduced. 

We have checked that the results (fT^ - (fT5|) and if^ de- 
rived here for d = 4 are consistent with the corresponding 
results of Ref. 4] for d < 4, so that (P) and ^ tend to 
(PI|l and (UHl) for d ^ 4, respectively. 

In the recent numerical study of the depinning transi- 
tion of driven interfaces at the upper critical dimension 
in the random-field Ising model (RFIM) the logarith- 
mic corrections to the force-velocity characteristics were 
chosen in the form 



[v(/)//]'/^^|ln/|. 



(25) 
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The best fit to numerical data was obtained with (j) = 
0.40 ± 0.09. Taking into account that the numerical de- 
termination of the logarithmic corrections is difficult, this 
value is in a fair agreement with our exact result — 2/9. 
The reason of the discrepancy might be due to the fact 
that simulations are not carried out in the as ymp totic 
regime, In / ^ —1. Indeed, the simulations in [23 were 
performed for | ln/| 1 4.5. 

In the remaining part of this Report we will discuss 
the contributions of the subdominant logarithms, which 
appear in the two- loop order of RG. Using the results of 
the work J§| we find that the two-loop correction to the 
disorder correlator at d = 4 has the form 

^2i{z) = Kl'-i^{Qaf)-^'^p2{z{-il/a)-^'^), (26) 

where p2{z) obeys the following differential equation 

{(1 - p{z)){p2{z) p^m + l/2[(p(z) - l)p'\z) 
+piz)]}" + zp'^iz) + 4p2iz) ^ 0, (27) 

with p{z) given by (|12|l . Following Q we impose the 
boundary condition /52(0) = {p2{z) is the counterpart 
of the function ?/2 {u) of work , so that this condition is 
required for the consistency of the results for both d < 4 
and d = 4). Expanding P2{z) in a Taylor series we obtain 
P2{z) — ~\z\ + 19/18z^ + .... The leading correction to 
Eq. (|T7| has the form -A;(0+)V(7g)'* Xt l/fc^(fc + q? 
0, so that the logarithmic correction to the roughness 
l)18|l behaves as (InxAo)"^/^ and, therefore, is irrelevant 
for large x. In order to take into account the two- loop 
corrections to the force-velocity characteristics we need 
the two-loop correction to the friction, which reads Q 

Klh^ (A"(0+)2 + A"'(0+)A'(0+)[3/2 - In 2]) . (28) 

Substituting Eq. (j^SJl into Eq. and Eq. I^T]) into 
Ea. (|28|l we obtain the two- loop contribution to Eq. H19|) 
as C/^^ where C = -[17/2-Kln2]/54 « -0.17. The two- 
loop correction to the driving force © results only in a 



shift of the threshold force Fc, which is a non- universal 
quantity, so that we will not consider it. Using the same 
arguments which led us to Eqs. H23|) and (|24|l we arrive 
at 

/~;y'exp(-2/v), v(/)//^/2/9(l + C/Zv). (29) 

Eqs. H29|l express the force- velocity characteristics to two- 
loop order of FRG by using the parametric representation 
in terms of the correlation length — InLvAg. The 
latter allows us to avoid the asymptotic expansion, which 
leads to a complicated expression [Tsj. 

In order to describe the crossover to the asymptotic be- 
havior given by Eq. (|24l) , one needs in addition to (|29|l the 
corrections to scaling. The latter are not available, since 
so far the only asymptotic solution of the flow equations 
©-Q is known. Due to these reasons the extension of 
the description of the critical behavior at the depinning 
transition to larger v is a non trivial problem. 

In conclusion, we have considered the effects of 
quenched random field disorder on the driven elastic in- 
terface at the upper critical dimension close to the de- 
pinning transition. We have shown that the interface 
undergoes the depinning transition at the critical driving 
force Fc, and we have obtained the logarithmic correc- 
tions to the displacement correlation function, the cor- 
relation length, and the force-velocity characteristics. In 
approaching the depinning transition our one-loop results 
become exact. We hope that the analytical results de- 
rived here will be of interest for numerical studies of the 
depinning transition. 

After this paper was finished, we learned that the more 
precise estimation of (j) from the simulations of RFIM 
is in a very good agreement with our prediction 2/9. 

A support from the Deutsche Forschungsgcmcinschaft 
(SFB 418) is gratefully acknowledged. We also would like 
to thank P. Le Doussal for a useful discussion. 



[1] D. S. Fisher, Phys. Rep. 301, 113 (1998). 

[2] H. Leschhorn, T. Nattermann, S. Stepanow, and L.-H. 

Tang, Ann. Phys. (Leipzig) 6, 1 (1997). 
[3] M. Kardar, Phys. Rep. 301, 85 (1998). 
[4] T. Nattermann, S. Stepanow, L.-H. Tang, and H. 

Leschhorn, J. Phys. 11 France 2, 1483 (1992). 
[5] O. Narayan and D. S. Fisher, Phys. Rev. B 48, 7030 

(1993). 

[6] P. Chauve, T. Giamarchi, and P. Le Doussal, Phys. Rev. 

B 62, 6241 (2000). 
[7] P. Chauve, P. Le Doussal, and K. J. Wiese, Phys. Rev. 

Lett. 86, 1785 (2001). 
[8] P. Le Doussal, K. J. Wiese, and P. Chauve, Phys. Rev. 

B 66, 174201 (2002). 
[9] H. Ji and M. O. Robbins, Phys. Rev. B 46, 14519 (1992). 
[10] U. Nowak and K. D. Usadel, Europhys. Lett. 44, 634 



(1998). 

[11] L. Roters et al, Phys. Rev. E 60, 5202 (1999). 

[12] A. Rosso, A. K. Hartmann, and W. Krauth, Phys. Rev. 

E 67, 021602 (2003). 
[13] H. Leschhorn, Physica A 195, 324 (1993). 
[14] M. Jost and K. D. Usadel, Physica A 255, 15 (1998). 
[15] A. Rosso and W. Krauth, Phys. Rev. Lett. 87, 187002 

(2001). 

[16] L. Roters, S. Liibeck, and K. D. Usadel, Phase Transi- 
tions 75, 257 (2002). 

[17] R. Chitra, T. Giamarchi, and P. Le Doussal, Phys. Rev. 
B 59, 4058 (1999). 

[18] J. Zinn- Justin, Quantum field theory and critical phe- 
nomena (Clarendon Press, Oxford, 1996). 

[19] N. N. Bogoliubov and D. V. Shirkov, Introduction to the 
Theory of Quantized Fields (Nauka, Moscow, 1984). 



L. Rotors, S. Liibcck, and K. D. Usadel, Phys. Rev. E [21] S. Liibeck (private communication). 
66, 026127 (2002). 



